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I,  Introduction 

This  is  a  siunmary  of  work  performed  under  Contract  No.  W28~099-ac-172 
awarded  "by  the  Watson  Laboratories  of  the  United  States  Air  Force  to  a  group  of 
mathematicians  associated  v/ith  the  Institute  for  Mathematics  gind  Mechanics  of  Nev/ 
York  University,  The  group  tegan  to  be  organized  in  July,  I9U6  to  conduct  basic 
mathematical  research  in  the  field  of  transmission,  reflection  and  diffraction  of 
high  frequency  electromagnetic  v/aves.  The  investigations  to  be  carried  out  were 
conceived  as  primarily  theoretical  in  nature.  However,  it  was  a  guiding  thought 
that  the  studies  be  advanced,  where  possible,  to  a  point  of  potential  engineering 
applicability, 

A  number  of  Research  Reports  were  submitted  by  the  group  to  the  Watson 
Laboratories,  They  are  discussed  in  this  report,  and  are  listed  by  title  and 
author  at  the  conclusion.  Reference  to  them,  in  the  report,  is  made  by  the  occurr- 
ence of  numerals  in  square  brackets. 

The  plan  of  this  report  is  as  follows.  The  next  section  will  deal  with 
work  on  propagation  in  stratified  media.  Part  III  v;ith  the  scattering  of  propaga- 
tion by  obstacles.  Part  IV  will  give  a  rather  detailed  account  of  work  in  progress, 
on  a  general  theory  of  the  relation  of  geometric  optics  to  electromagnetism.   Inso- 
far as  this  v/ork  is  in  part  directed  to  the  discovery  of  asymptotic  expressions  for 
certain  integrals  it  represents  a  continuation  of  the  work  described  in  II  where 
such  integrals  are  used  (see  item  3  of  II ).  Hov;ever,  the  approximate  evaluation 
of  these  complex  integrals  is  only  a  connecting  link  betv/een  two  distinct  problems. 
In  V  there  are  brief  indications  of  other  work  of  the  group,  some  of  which  came  to 
be  abandoned  at  least  for  the  time  being. 

II,  Propagation  in  Stratified  Media. 

One  of  the  principal  problems  studied  by  members  of  the  group  was  the 
propagation  of  electromagnetic  radiation  in  a  space  with  variable  £ ,  yu  ,  and  O"  , 
This  problem  becomes  approachable  only  if  it  is  restricted  by  requiring  that  the 
electroraeignetic  parameters  be  functions  of  a  single  space  coordinate.  In  the  case 
of  a  plane-stratified  medium,  for  example,  one  assumes  that  £  ,  a<  ,  and  CT  vary 
with  one  coordinate  only,  and  the  medium  is  characterized  by  the  assignment  of 
three  functions,  f  (z),  m(z),  <5'(z),  where  z  is  one  of  the  rectangular  space- 
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coordinates,  A  problem  is  now  defined  "by  the  specification  of  a  source  of  electro- 
magnetic radiation,  e.g.,  a  harmonic  dipole  located  at  the  origin  of  coordinates. 
One  then  seeks  formulae  which  describe  the  radiation  field  everywhere  in  space. 
This  class  of  problems  embraces  the  so-called  problem  of  "anomalous  propagation": 
i.e.,  the  transmission  of  radar  intelligence  through  the  atmosphere  (account 
being  taken,  in  some  manner,  of  the  sphericity  of  the  earth). 

This  problem  has  received  a  great  deal  of  attention,  particularly  in 
recent  years  when  its  importance  in  radar  communication  became  evident,  and  many 
references  to  it  \/ill  be  found  in  the  Bibliography  |9j  prepared  by  members  of  this 
group.   In  spite  of  the  brilliant  work  already  devoted  to  this  problem,  the  great 
importance  of  the  problem  from  a  theoretical  as  v;ell  as  a  practical  point  of  view 
calls  for  further  effort,  particularly  in  view  of  the  extremely  formidable  numeri- 
cal calculations  to  which  this  work  gives  rise.  Moreover,  the  present  solutions 
utilize  approximations  which  seriously  limit  the  range  of  applicability  of  the 
results, 

1,   In  the  first  Report   j 3  |  submitted  on  this  general  subject,  it  is  assumed 
that  the  half-spaces  z  <  0  and  z  >X (a  positive  constant)  are  occupied  by  homogene- 
ous substances,  and  that  in  the  remaining  space,   0  <  z  <  /  ,  there  are  given  three 
functions  £(z),  /w(z),  O  {z)   which  are  piecewise  different! able  in  the  interval 
0  to  X  ,   It  is  assumed  that  a  plane  wave  approaches  the  inhcanogeneous  slab  from 
below,  and  the  problem  is  to  investigate  the  effect  vrhich  the  layer  system  has  on 
the  propagation  of  the  wave. 

This  problem  is  approached  through  the  special  case,  itself  of  great 
interest,  of  a  finite  multilayer.  Here,  the  interval  0  to  X  is  supposed  divided 
into  a  finite  nuraber  of  intervals:  0  =  jl      <  j(     <   i/     <  , , .    <  ^     =  /^  in  each  of  which 

^(z)f  /^(z),  <^(z)  are  constant,   For  this  case  there  is  derived  in  the  paper  a 
complete  set  of  recursion  formulas  (  [3]  ,  p.l7,  formulas  25  and  26)  permitting 
the  ready  calculation  of  all  desired  field  quantities.  The  case  that  £  ,  /u  ,  (T 
are  piece\/ise  continuously  different iable,  is  then  treated  by  a  passage  to  a  limit. 
Here  the  reflection  and  transmission  coefficients  of  the  layer  system  are  expressed 
relatively  simply  in  terms  of  certain  fanctions,  which,  however  are  determined  by 


Blccatl  equf-tions  vfhose  variable  coefficients  depend  on  the  parameter  functions 
£(z),  yA(z)  and  <5(z). 

2,  A  second  Report  [_Uj  considerably  generalized  the  multilayer  case  men- 
tioned above  by  introducing  a  dipole  source  of  radiation.  The  first  part  of  the 
paper  investigates  the  case  of  a  multilayer  contained  between  tv/o  semi-infinite 
homogeneous  media  v/ith  a  stream  of  plane  waves  of  constant  amplitude  impinging  on 
the  layer  from  below,  the  \/aves  being  so  polarized  that  the  magnetic  vector  is 
parallel  to  the  multilayer,  IText  a  dipole  source  below  the  multilayer  is  treated, 
the  source  being  represented  as  a  superposition  of  plane  waves.  Finally,  the 
dipole  is  located  v/ithin  the  multilayer.  The  special  case  of  a  single  layer  is 
worked  out  in  detail.   The  radiation  fields  are  expressed  in  the  form  of  contour 
integrals,  which  take  account  of  the  representation  of  the  dipole  source  as  a  super- 
position of  plane  v^aves. 

The  methods  of  the  paper  raise  certain  difficult  questions  of  a  mathe- 
matical nature  concerning  the  convergence  of  the  integral  representations.  This 
question,  as  v/ell  as  the  problem  of  effectively  calculating  v/ith  the  form  of  solu- 
tion given  in  the  p^er  were  postponed  in  favor  of  the  more  general  investigation 
described  belovr. 

3.  In  a  fundamental  paper  I  6  J  the  propagation  of  electrcaaagnetic  v/aves 
through  a  quite  general  plane  stratified  inhomogeneous  meditim  was  studied  for  the 
case  of  an  arbitrary  source.  It  v/as  shown  that  the  problem  may  be  reduced  to  the 
study  of  a  pair  of  scalar  functions,  each  subject  to  a  partial  differential  equa- 
tion, the  two  equations  differing  only  in  the  interchange  of.  £.'(z)  and  /^(z,),  where 
these  occur  (  }_6jp.l2,  formula  U.7).   For  fields  periodic  in  time,  the  problem 

is  thOn  reduced  to  the  superposition  of  solutions  of  differential  equations  of 
Riccati  type  t  [_6jp.  J>2   formulas  (7.IO)  and  (7. 11))  .,  This  formxilation  has  the  ad- 
vantage that  the  solution  functions  of  the  Riccati  equations  express  very  directly 
the  distortion  in  \/ave  front  and  the  change  in  amplitude  of  a  plane  wave  incident 
on  a  stratified  medium.  The  total  electromagnetic  field  in  the  medi\im  is  expressed 
as  a  complex  integral  involving  the  representation  of  the  arbitrary  source  as  a 
superposition  of  plane  v/aves,  and  the  solution  functions  of  the  Riccati  equation 
(^[e]   .formulas  (7.8)  and  (7.9)}  » 
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The  formulas  of  this  paper  constitute  a  solution  of  the  general  prohlem 
in  the  sense  that  it  is  possible  to  apply  to  them  numerical  procedures  to  calculate 
the  fields  associated  v/ith  a  specified  source  and  chosen  set  of  parameter  functions 

£(z),  //*(z),  ^(z).  However,  some  initial  efforts  in  this  direction  indicate 
that  the  task  would  "be  an  extremely  formidahle  one,  and  there  is  no  reason  to 
■believe,  as  yet,  that  such  an  undertaking  vroiild  offer  advantages  over  the  current 
procedure  of  "normal-mode"  calculations  based  on  methods  of  Furry,  Pekeris,  and 
others.  (See,  among  many  other  references  listed  in  the  Bibliography  [9],  The 
Radiation  Laboratory  Report  ^=680,  Theory  of  Characteristic  Functions  in  Problems 
of  Anomalous  Propagation,  Februaiy  28,  19^+5) • 

It  is  felt  that  nuich  further  work  of  a  theoretical  nature  must  first  be 
carried  through  in  connection  with  approximate  evaluations  of  the  complex  integrals 
appearing  in  this  paper,  and  it  is  hoped  that  the  project  described  in  Part  IV  of 
this  report  may  supply  some  of  the  methods  for  such  an  evaluation* 

k.  Since  the  initial  motivation  for  the  studies  described  above  v/as  the 

problem  of  "anomalous  propagation",   some  thought  has  been  given  to  the  problem  of 
adapting  a  theory  for  a  plane-stratified  medium  to  the  case  of  the  earth.  This  is 
conventionally  handled  by  the  use  of  a  modified  index  of  refraction.  The  problem 
is  discussed  in  the  Introduction  to  [Jo]   ,  reported  on  above,  and  in  [8j  ,   It  will 
be  shown  in  a  later  Report  that  the  same  modified  index  of  refraction  can  be 
applied  to  the  "flat-earth"  theory  of  [6j  .  However,  this  device  is  not  satis- 
factory in  some  applications,  and  it  is  desirable  that  the  problem  be  approached 
in  some  other  \;ay.   The  most  direct  approach  would  seem  to  be  an  investigation  6-^ 
spherically  stratified  media. 

One  such  study  is  begun  in  [sj  v/hich  deals  with  a  spherically  symmetric 
medium.   Here  the  electromagnetic  parameters  depend  only  on  the  distance  to  the 
origin  of  space.  This  study  parallels  the  first  part  of  a  preceding  Report  [s]   by 
showing  that  the  components  of  an  electromagnetic  field  in  a  mediiim  v/hose  structural 
constants  depend  only  on  the  distance  to  the  origin  can  be  derived  from  two  scalar 
functions.  The  use  of  a  special  curvilinear  coordinate  system  simplifies  the  form 
of  the  differential  equations  which  govern  these  scalar  functions  (£8^  p.l6,  (5.9^ 
and  makes  easier  the  calculation  of  the  components  of  the  electromagnetic  field 
(ibid, (5.10 )j  .  However,  at  the  point  at  which  this  Report  stops,  the  differential 
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equations  are  not  explicitly  solved.   The  Report  brings  out  the  striking  similari- 
ties "between  the  equations  associated  with  the  spherically  stratified  and  the  plane 
stratified  media.   It  is  hoped  that  some  advantage  can  he  taken  of  this. 

5.  Another  approach  to  the  problem  of  anomalous  propagation  over  a  spherical 
earth  is  no\;  \/ell  under  way  carrying  over  to  the  case  of  spherical  stratification 
the  method  of  characteristic  functions  (see  the  previously  cited  Radiation  Labora- 
tory Report)  ncM   employed  for  a  flat  earth  and  horizontally  stratified  atmosphere. 
The  object  of  this  investigation  v;ill  be  to  determine  v/hether  this  more  direct 
^proach  to  the  problem  of  "anomalous  propagation"  in  the  earth's  atmosphere  v/ill 
substantially  complicate  the  task  of  numerical  calculation  associated  with  the 
"flat-earth"  theory. 

6.  Reference  has  already  been  made  to  work  in  progress  on  the  problem  of 
the  numerical  calculation  of  anomalous  propagation  fields  on  the  basis  of  the 

formalism  developed  in  [6  j  .  The  work  is  being  undertaken  by  one  of  the  members 
of  the  group  well  acqiiainted  with  calculation  procedures  used  in  conjunction  v/ith 

modern  computing  instrumentation.  The  object  of  the  study  is  to  determine  whether 

direct  computation  of  the  integrals  found  in  [^6jseems  practically  feasible.  Results 

so  far  are  inconclusive,  and  the  investigation  is  being  continued, 

III.  The  Scattering  of  Radiation  by  an  Obstacle. 

The  general  problem  of  scattering  of  radiation,  i.e.,  the  reflection 
from,  transmission  through  and  diffraction  by  an  arbitrary  obstacle,  can  be  sub- 
sumed under  the  larger  topic  of  propagation  in  non-homogeneous  media,-  However, 
the  problem  is  here  construed  in  uhe  narrower  sense  that  an  obstacle  of  one  kind 
of  material  is  imbedded  in  an  otherv/ise  homogeneous,  electrically  isotropic  space 
in  v/hich  there  is  present  some  source  of  electromagnetic  radiation,  far  removed 
from  the  source..  The  condition  of  large  distance  of  the  source  from  the  obstacle 
is  made  to  insure  that  the  field  reflected  back  to  the  source  is  negligible  in 
comparison  uith   the  field  initiated  there.  The  problem,  given  the  nature  of  the 
source,  the  aedi\im,  and  the  obstacle,  and  the  geometry  of  the  configuration;  source 
and  obstacle,  is  to  find  the  total  field  at  all  points  of  space..  This  problem  can 
be  solved  exactly  only  in  cases  where  the  obstacle  is  of  particularly  simple 
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geometric  character.  In  recent  years  several  papers  have  appeared  which  have  in- 
vestigated the  problem  for  the  case  of  fairly  arbitrary  obstacles,  using  approxi- 
mate methods  (see  references  in[_5j,  p.  ^)-  Approximate  methods  of  a  very  general 
application  to  problems  of  this  kind  consist  in  the  methods  of  geometrical  optics, 
corresponding  to  the  limit  of  zero  v/ave  length,  and  the  Kirchhoff  method,  which 
will  be  briefly  indicated  below.  Since  these  methods  rest  upon  approximations  to 
electromagnetic  theory,  results  obtained  by  them  ought  to  be  checked  against 
methods  which  are  more  accurate  in  the  boundary  conditions  imposed  upon  the  fields 
at  the  surface  of  the  obstacle.   Such  an  investigation  is  undertaken  in|_5j,  dis- 
cussed below, 

1,    However,  in  the  first  Report [^ij on  this  general  problem  the  method  used 
is  the  Kirchhoff  method,  and  some  of  the  results  obtained  by  it  are  checked  sigainst 
the  method  of  geometrical  optics.  The  obstacle  is  assumed  to  be  perfectly  reflect- 
ing, but  of  arbitrary  smooth  shape.   The  wave  equation  satisfied  by  the  total 
vector  field  is  transformed  to  an  integro-differential  equation.   In  general,  this 
eqiiation  involves  both  surface  and  volume  integrals.   In  the  present  case  of  a 
perfectly  reflecting  obstacle,  the  volume  integral  is  absent.  Nov;  the  essential 
featiire  of  the  Kirchhoff  method  is  invoked.  This  is  the  assumption  of  boundary 
conditions  which  are  false  in  the  electromagnetic  case,  but  piausible  in  geom- 
etrical optics:  e.g.,  it  is  assTimed  that  the  total  field  is  double  the  incident 
field  on  the  geometrically  illvuainated  part  of  the  surface,  and  zero  elsewhere 
on  the  surface.  The  effect  of  these  assumptions  is  to  reduce  the  problem  of  finding 
the  total  field  to  the  problem  of  evaluating  a  definite  surface  integral.  This 
evaluation,  while  theoretically  possible  given  the  geometric  shape  of  the  obstacle, 
is  prohibitively  difficult  to  perform  exactly.  An  approximate  evaluation  is  found, 
essentially  by  the  method  of  stationary  phase  applied  to  surfaces. 

Effectively,  the  succeeding  Report j 5] is  of  greater  generality,  and  a 
discussion  of  the  nature  of  the  results  obtained  herelljwill  be  postponed  to  the 
next  section. 
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2,    It  is  possible  to  set  up  a  pair  of  integro-differenti?!  eqaations  for 
the  electric  and  magnetic  fields  in  the  presence  of  a  general  ohstacle,  without 
the  restriction  that  the  obstacle  te  perfectly  reflecting.  But  without  sone 
restriction  on  the  obstacle,  these  equations  are  completely  unmanageable.   It  was 
observed  in  the  acoustic  case  by  H.  Prinakoff  and  J.  B.  Keller  (see  ref.  12  in 
|5 J  )  and  carried  over  to  the  electromagaetic  ease  inr5j  that  these  equations  can 
be  reduced  for  smooth  curved  obstacles  of  constant  thickness  h,  if  h  is  sufficient- 
ly small.   In  this  case,  the  volume  integrals  may  be  replaced  by  surface  integrals, 
and  certain  differentiations  of  field  quantities,  can  be  carried  out  e3q)licitly. 

The  boundary  conditions  imposed  are  those  of  electromagnetism,  so  that 
the  approximations  made  in  geometrical  optics  or  in  the  Kirchhoff  method  are  now 
avoided.  There  remains  the  approximation,  as  in  [l  [that  a  certain  surface  integral 
must  be  evaluated  by  some  approximate  method.  The  integral  is  of  the  general  form: 

f  ,   i(ki|r'-r  I  t-  k  |r'-r  J  ) 

J^     f(r«)  e   "■     ^     "^     '^     ds 

\fhere  r'  is  the  vector  position  of  a  variable  point  on  the  surface  S,  r  and  r 

locate  the  source  point  and  observation  point,  f(r')  is  a  certain  function  slowly 
varying  over  the  surface  in  comparison  with  the  exponential,  and  k-,  and  k  are 

certain  constants.   This  can  be  evaluated  by  an  adaptation  of  the  method  of 
stationary  phase. 

The  solution  fields  are  obtained  in  the  foru  of  an  expansion  in  po'ers 
of  h  and  the  first  two  terms  can  be  explicitly  calculated  by  the  methods  of  this 
paper.   The  first  term  is  merely  the  incident  field,  i,e.  the  field  due  to  the 
source  in  the  absence  of  any  obstacle.  The  first  correction  term  to  the  incident 
field  is  sho\m  to  e^foal  the  incident  field  of  the  source  multiplied  by  three  fac- 
tors, designated  as  the  geometric,  the  phase,  and  the  reflectivity  factors.  The 
first  erf  these,  in  the  limit  of  zero  wave  length  agrees  v/ith  the  reflection  co- 
efficient obtained  by  geometrical  optics,  and  the  second  with  a  corresponding 
phase  factor  which  results  from  an  application  of  the  Kirchhoff  method.  The 
third  factor  is  not  present  in  the  JCirchhoff  solution  and  represents  an  actual 
improvement  in  the  accuracy  of  the  solution. 
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The  paper  [5 J  discusses  the  physical  significance  of  these  factors  in 
considerahle  detail  and  conpares  then  v/ith  results  previously  obtained  ty  other 
methods,  (see,  in  particular55  »  ihid, ,  pp.  26-35)«  The  agreement,  in  all  cases, 
appears  to  be  satisfactory.  As  a  further  check  upon  the  theory  the  explicit  for- 
inalas  were  derived  for  the  exact  fields  corresponding  to  the  geometrically  simple 
case  of  a  spherical  shell  \/ith  a  radiatin^;  dipole  at  the  center.  This  is  dis- 
cussed below, 

3»     In  a  subsequent  Report|_7jexplicit  expressions  are  obtained  for  the 
exact  fields  every\/here  in  space  associated  with  a  radiating  dipole  at  the  center 
of  a  sphere  of  one  nedium  encased  in  a  spherical  shell  of  raediun  2  \/hich  is  bounded 
by  mediiiEi  3«  The  electromagnetic  parameters  are  constant  within  each  medium,  but 
othen/i&e''arbitrary,  and  the  shell  radii  are  arbitrary.  The  reflected  and  trans- 
mitted fields  are  examined  in  various  special  cases  and  a  variety  of  results  knovm 
for  flat  media  are  obtained  by  limiting  considerations  applied  to  the  shell  radii. 
Comparison  is  made  of  the  exact  results  of  this  special  case  with  the 
approximations  obtained  by  the  general  theory  above  \'5j*     I't  is  found  that  except 
for  one  field  component  and  except  near  the  singularity  of  the  reflected  field,  the 
approximate  solution  is  good, 

U.     A  study  is  under  way  of  the  general  problem  of  effectively  calculating 
by  the  formulas  of  [5j  the  fields  associated  with  a  particular  source  and  obstacle. 
Since  the  formulas  in  \^j   are  quite  complicated,  necessairily  so  as  one  sees  by  con- 
sidering the  complexity  of  the  exact  fields  associated  with  the  simple  special  case 
investigated  in  Uj   ,  it  is  natural  to  inquire  hov;  usef;il  these  formulas  may  be 
for  quantitative  results.   The  study  is  confined  to  the  case  of  quadric  surfaces  and 
special  positions  of  the  radiating  source. 

5,  Another  study  in  progress  seeks  to  develop  formulas  for  the  case  of  an 
interface  of  arbitrary  shape  between  t\/o  semi-infinite  media..  This  work  was  in 
abeyance  in  the  past  months;   it  is  described  in  some  detail  in  the  preceding 
Q;aarterly  Progress  Report, 

6.  A  quite  different  approach  to  scattering  problems  is  nearing  completion, 
using  variational  methods.  This  project  is  devoted  to  calculating  the  back- 
scattering  cross-section,  that  is,  the  power  radiated  back  v;hen  a  plane  v,rave  of  .unit 
power  density  per  unit  area  impinges  on  a  target,  for  a  variety  of  obstacles. 


9. 


The  differential  equations  and  exact  "boundary  conditions  of  the  iDro'blem 
are  converted  to  integral  equations.   It  is  shoun  how  the  protlem  of  solving  these 
equations  may  he  transformed  into  the  prohlem  of  finding  extreme  values  of  certain 
integrals.   This  prohlem  is  treated  by  the  methods  of  the  calculus  of  variations. 

The  first  prohlem  treated  in  this  way  is  the  scattering  of  sound  hy 
rigid  bodies  of  arbitrary  shape.  Here  the  back-scattering  is  found  directly  as 
the  extreme  value  of  the  ratio  of  two  integrals.  Examples  are  calculated  for  cer- 
tain simple  shapes  v;here  exact  answers  can  be  obtained  in  order  to  estimate  the 
accuracy  of  the  method  of  the  paper.   Another  problem  treated  is  the  scattering  of 
electric  waves  by  cylindrical  obstacles  where  the  cross  section  of  the  cylinder  may 
be  arbitrary.   Finally,  it  is  indicated  \;hat  other  problems  may  be  treated  by  the 
methods  of  the  paper. 

To  bring  out  the  details  of  the  method  one  may  consider  the  scattering 
of  sound  by  obstacles  of  arbitrary  shape.  Denote  the  velocity  potential  by 

— ikct 
V^(x,y,z)  e      ,  which  satisfies  the  v/ave  equation: 

1)   (^^  +  k^)  y^(x,y,z)  =  0 
The  obstacle  is  taken  as  rigid  and  so  the  normal  derivative  of  if 
vanishes  on  it.  Assume  that  the  field  incident  on  the  obstacle  may  be  represented 
as  an  incoming  plane  v/ave,  propagated  in  the  positive  x-direction.  The  total  field 
contains,  also,  an  outgoing  spherical  wave,  whose  angular  pattern  is  to  be  deter- 
mined. Let  the  value  of  ^  on  the  obstacle  be  denoted  by  *^  (a)*.  Then  it  is  shown 
that  .^  satisfies  the  integral  equation: 

^'g,  Ir-r'l 

Here  the  integral  extends  over  the  surface  of  the  obstacle,  and 

Ir-r'i   =  (x-x')  +  {-^-Y^)     +  (z-z')  ^     The  symbol  — ^  denotes  the  normal 

5  n^ 

derivative.  The  condition  that  the  normal  derivative  of  ^vanishes  on  the  obstacle 
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gives  the  integral  equation,  valid  on  the   obstacle: 


3) 


a      ikx 

0  =  — —     e 

an 


.-  -ferLv'^'^ST 


9      e 


.kjr-r'l 


^n 


Ir-r'l 


ds' 


Here  S  and  S'    refer   to   the   surface   of   the   ohstacle,    the  primed  co- 
ordinates  refer  to  a  running  point   on  the   obstacle,    the  unprimed  tea  fixed  point, 
This   is   an  integral   equation  of   the  form: 


U) 


H(s)  =    /      <^(s')   K(s,s')   ds« 


The  general  theorem  can  be  proved  that  the  function  v/hich  renders 
stationary  the  integral  expression: 


5)  Ktf  )  = 


,^(s)  K(s,s')</?(s')  ds  ds' 


[|(/?(s)  H(s)  dsl 


is  the  function  which  satisfies  the   integral   equation  U)   save   for  a  multiplicative 
constant  v/hich  may  he  used  to  normalize  ^to  unit    incident  pov/er  density.      Con- 
sequently,   the   equation  3)  ^^y   '^e   expressed  in  the  variational    form   5). 

If  \/e  compute  the  asymptotic  form  of  tp  from  the  integral  expression  2) 
we  obtain: 


ikx 


ikr 


6)  Cf  -^  e^'^  +  f{Q,(f) 


Here  r,  © ,  (^   are  spherical   coordinates  corresponding  to  x,y,z.      The 
function  f(  6  ,  W)  becomes: 

7)         f(0,(^)  =  _     1       /(^.(st)  _i_  g-ikr  •  [cos  0  cos  e  i+sinOsin  ©  icos(  ^'~  V?' )]  ^^  , 
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Here  r'   0'    ^'   are  s-oherical   coordinates   on  the  surface   S'.      One  may 

'  2 


interpret  f{9, Cf)   "by  saying  that  dil  f{Q,i/>)\       is  the  power  scattered  into  an 
element  of  solid  angle  dO^  around  the  direction  9,  (J^  .  The  "backvmrd  direction  is 
given  "by  0  =  TT,  for  v;hich  f(rr,  (^  )   is  independent  of  t^  and  equal  to: 

This  expression  may  be  shovm  to  tJt?  equal  v;ithin  a  constant  to  l((^), 
so  the  back-scattered  amplitude  may  be  calculated  directly  by  considering  l(^) 
stationary.  This  may  be  done  by  standard  methods  of  the  calculus  of  variations, 

IV,  Asymptotic  Developments  of  Steady  State  Electromagnetic  Fields. 

1,     During  this  past  year  new  investigations  uere  started  on  the  subject 
of  asymptotic  expansions  of  steady  state  electromagnetic  fields.   These  expansions 
express  the  steady  state  field  in  the  form  of  a  series,  not  necessarily  convergent 
but  yielding  the  field  fairly  closely  if  even  only  two  or  three  terns  of  the  expan- 
sion are  used,  provided  however  that  the  basic  variable  in  the  expansion  be  large. 
The  basic  variable  in  the  expansions  sought  in  this  research  is  the  reciprocal  of 
the  \/avelength,  so  that  for  small  v/avelengths  the  first  few  terms  of  the  expansion 
might  give  a  close  approximation  to  the  steady  state  field. 

Asymptotic  expansions  are  not  nev/  to  electromagnetics  but  their  applica- 
tion is  thus  far  not  extensive.   It  is  also  not  a  new  idea  to  use  geometric  optics 
solutions  as  approximations  to  steady  state  electromagnetic  fields  for  small  wave 
lengths,  vrhat  is  nev;  in  this  investigation  is  the  development  of  a  comprehensive 
theory  which  v/ill  furnish  asymptotic  expansions  in  which  the  geometric  optics  solu- 
tions appear  as  the  first  term  of  the  expansion.  This  theory  v;ould  then  show  the 
validity  of  the  geometric  optics  solution  in  any  particular  application,  and  by  the 
use  of  more  terns  of  the  series  afford  a  better  approximation  to  the  true  field. 

The  follov/ing  result  is  basic  in  this  investigation.  Consider  a  source 
of  electromagnetic  waves  which  has  any  distribution  in  space  but  which  at  time 
t  =  0  suddenly  rises  from  0  amplitude  to  amplitude  1,   It  is  not  supposed  that  space 
is  homogeneous.  The  field  satisfying  Maxwell's  equation  and  arising  from  this  source 
is  called  the  pulse  solution  of  Maxv/ell's  equations  or  the  pulse  field.  If  one  now 

fixes  on  a  point  in  space  there  will  be  discontinuities  in  this  pulse  field  at  this 
point  at  particular  tines  t^ .  These  are  the  instants  at  which  either  the  directly 
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transmitted  wave  from  the  source  or  v/aves  reflected  from  obstacles  in  space  pass 
this  point,.  The  discontinuities  arise,  of  course,  from  the  fact  that  at  any  time 
t^y  a  J^6v;  wave  passes  this  point  v/hose  amplitude  is  zero  at  time  less  tlian  t^^ 
hut  some  non-zero  value  at  times  greater  than  t^.^.   There  \/ill  also  be  discontinu- 
ities in  the  derivatives  of  the  pulse  solution  as  well  as  in  ttie  pulse  solution  it- 
self. The  hasic  result  is  tliat  the  first  term  in  the  asyiaptotic  expansion  of  any- 
steady  state  electromagnetic  field  arising  from  a  harmonic  source  having  the  same 
space  distrihution  as  the  pulse  is  given  hy  the  discontinuities  in  the  pulse  solu- 
tion itself.  Moreover,  the  successive  terms  in  the  asymptotic  expansion  are  given 
hy  the  discontinuities  of  the  successive  derivatives  of  the  pulse  solution.  Ex- 
plicitly these  formulas  are  as  follov/s: 

00 


u=     -— g 


dE.  _  ,,,.,;       ,        r  £^ 


[El     e'-'Y.- 4^11   [4^1      e"=^--^Al4e'^''dt 
L   qJcx  X  u)  ^ —    '-at    -'^  xuij      ^     2 


00 


.^ 


L    oi^  lOJ  ^—   L  at    J^  x(^)J      ^^ 

In  these  formulas,  u  and  v  are  the  vector  electric  and  magnetic  field 
respectively,  g  is  the  space  distribution  of  an  inqjressed  electric  field  of  the 
form  g(x,y,z)f(t),  and  E.  and  H  are  the  electric,  respectively  magnetic  fields  of 
the  pulse  solution.   The  square  brackets  with  a  subsceipt  o<  denote  the  jTimp  in  the 
enclosed  quantity  at  time  t_^.   The  functions  ^y/^give  the  v;ave  fronts  at  each  point 
of  space  at  successive  times  t ^through  the  condition  that  the  equations 
ct  =  yjQ^iXfYfZ)   are  the  hypersurfaces  in  space-tine  corresponding  to  \/ave  fronts. 

The  connection  between  the  discontinuities  in  the  pulse  solution  and 
geometrical  optics  is  essentially  that  the  discontinuities  are  those  furnished  by 
v/ave  fronts  and  the  behavior  of  the  wave  ^  ront  of  aiiy  signal  is  the  behavior  of  the 
electromagnetic  field  as  furnished  by  geometrical  optics.  This  connection  which 
is  here  utilized  is  not  new  (see  Luneburg,  E, :  Mathematical  Theory  of  Optics, 
Brovm  University,  I9UU)  .  In  the  paper  nov/  being  planned  on  the  theory  of 
asymptotic  solutions  some  of  this  older  material  v/ill  be  included  for  completeness 
sake  and  for  ready  reference. 


13. 


2,     The  following  material  represents  a  substantial  sketch  of  a  major  appli- 
cation that  v/ill  "be  made  of  the  theory  just  described;   that  is,  it  v/ill  be  shown 
how  the  asymptotic  field  transmitted  by  a  lens  can  be  obtained.  The  vord  lens 
here  can  be  interpreted  in  the  usual  optical  sense.   The  theory,  hov/ever,  is  broader 
in  that  it  applies  to  ani''  medium  bounded  by  tv/o  parallel  planes  separating  tv/o  semi- 
infinite  homogeneous  media. 

Since  the  general  theory  relates  the  asymptotic  behavior  of  a  steady 
state  electromagnetic  field  to  the  discontinuities  of  the  pulse  solution  and  the 
first  term  of  the  development  is  essentially  the  geometrical  optics  solution  to  the 
field,  it  is  \inderstandable  that  the  theory  finds  ready  ^plication  to  problems  in- 
tensively studied  in  the  well  developed  field  of  geometrical  optics,  v/here  this  much 
of  the  solution  to  the  steady  state  field  is  laiown,  A  brief  sketch  of  the  applica- 
tion nov/  follows; 

It  is  to  be  understood  that  the  material  in  the  sketch  presented  below, 
while  fairly  extensive  in  itself,  still  requires  considerable  development,  and  that 
a  final  paper  on  the  subject  may  not  be  forthcoming  in  quite  a  few  months. 

2.1  Let  us  first  consider  the  problem  of  transmission  of  light  through 
optical  instruments  from  the  point  of  view  of  the  electromagnetic  theory  of  light. 
The  x,y,z  space  is  divided  into  two  half  spaces,  called  object  and  image  space. 
These  two  half spaces  z  <  z  and  z  >  b  are  separated  by  the  optical  instrument , i.e. 

by  a  medium  z  <  z  <  z^  in  v/hich  the  material  constants  £  ,  M  ,  (5"  vaiy  with  x,y,z. 

We  assume  that  £  =  /^=  1,0"=  0  in  object  and  image  space. 

A  time  periodic  point  source  is 
located  at  a  point  P  of  the  object 
space.   The  problem  is  to  find  the 
P"°^^\;^y^/-;^^        ^T"*"^"       time  periodic  electromagnetic  field 


Pig.  1 


S  =  u(x,y,e)  e 
H  =  v(x,y,z)  e 


(2.11) 
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which  is  established,  in  the  image  space,  under  the  influence  of  this  source.  The 
vectors  u  and  v  then  are  solutions  of  the  time  free  Maxwell ^s  equations 


#.  -^ 


curl  V  +  iku  =0  ,    . 

^  ^  (2.12) 

curl  u  -  ikv  =0 

i»)  _  2n 

wnere  k  = r — 

c    A 

It  is  clear  tiiat  the  solution  u,v  of  our  problem  must  be  somewhat 
related  to  the  geometrical  optics  of  our  instrument,  in  particular  to  the  system 
of  light  rays  which  originates  at  the  point  source  P.  This  relation  must  be  such 
that  for  small  wave  lengths,  X,    i.e.,  for  large  values  of  k  =  r—  ,  the  vectors  u,v 

A 

shov;  an  asymptotic  behavior  v;hich  expresses  the  geometric  optical  situation  pro- 
vided by  the  instrument.  In  other  v/ords  we  are  concerned  vrith  a  certain  problem 
of  correspondence:  To  find  a  solution  u,v  of  (2.12)  which,  for  k  -»  oo  ,  shov/s  a 
prescribed  asymptotical  behavior, 

2.2  We  characterize  the  rays  of  the  instrument  by  Hamilton's  mixed 
characteristic  WCx  ,y  ,z  ;  p,  q).  This  is  a  function  of  the  coordinates  x  .y  .z 

— ■ 0   0   0  0   0   o 

of  the  source  and  of  the  direction  cosines  p,q,r  of  the  rays  in  the  image  space. 
It  determines  the  optical  length  of  the  ray  between  the  point  P     and  the  base  point 

Q,  of  the  perpendicular  dropped  from  the   origin  0  onto  the  ray  p,q,r.     We  can  derive 
all  essential  characteristics  of  the  ray  system  v/ith  the  aid  of  V.     The  ray  of 
direction  p,    q,   r,    for  exanple,    is  given  by  the  equation 

x-Jz+W       =0 

^  (2.21) 

y  ~  ^  z     +     U       =0 

If  we  calculate  p  and  q  from  these  equations  as  functions  of  x^y,z  and 
introdu.ce  the  result  into  the  expression 


^  =  px  +  qy  +  rz  +  W  (2.22) 
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then  a  function  ^=  ^  (x,y,z)  is  obtained  v/hich  gives  the  optical  length  of  a  ray 
connecting  tne  source  P  \/ith  the  point  x,y,z.   The  surfaces  y (x,y,z)  =  constant 

thus  are  the  \/ave  fronts  of  the  system  i.e.,  surfaces  orthogonal  to  the  rays. 

There  may,  of  course,  exist  several  v;ave  fronts  through  the  same  point 
x,y,z  just  as  there  may  exist  several  rays  \/hich  connect  this  point  with  the  source 

Let  dM  he  the  surface  element  in  v;hich  a  narrow  "tube"  of  rays  inter- 
sects a  wave  front  Y=  constant.   We  measure  the  expansion  of  the  ray  system  by 
the  expression 

^=  '^^^1_2_2  (2.23) 

which  in  our  case  of  a  homogeneous  space  is  equal  to  the  Gaussian  curvature  of  the 
wave  front. 

In  geometrical  optics,  light 
rays  are  considered  as  the 
carriers  of  light  impulses. 
If  the  light  distribution  along 
a  ray  is  observed  then,  accord- 
Pig,  2  ding' to  geometrical  optics,  the 

polarisation  of  light  is  unchanged  whereas  intensity  and  energy  flux  vary  propor- 
tional to  the  above  quantity  K.  We  may  describe  these  observations  mathematically 

by  a  vector  f ield  U  (x,y,z);   the  direction  of  these  vectors  determines  the 

I  |2 
polarisation  of  the  light,  the  absolute  value  |U|  the  energy  and  flux.  The 

vectors  u   along  a  fixed  ray  p,  q,  r  are  normal  to  the  ray  and  are  related  by  an 

equation  ->   -^ 
II  U  =  A(p,q)s^     (2.2U) 

where  the  vector  A(p,q)  is  constant 
along  the  ray. 

Fig.  3 
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2,3     The   actual   solution  of   the  above   transmission  problem  reauires 
the   integration  of  Maxwell's   equations 

curl  V     +  ik  S.  u     =0 

^  ^  (2.31) 

curl  u     "  ik  /^  V     =0 

v;here  £  and  u  are  given  functions  of  x,  y,  z.  On  account  of  the  npthematical 
difficulties  which  oppose  such  an  approach  one  prefers  to  establish  solutions  in 
the  image  space  (  S  =^=  l)  by  a  hi'pothesis: 

1)  The  vectors  u  and  v  are  analytical  solutions  of  (2,12)  which  are  regular 
and  uniformly  bounded  in  the  total  x,y,z  sijace, 

u       — ik  U^  V       — ik  ^ 

2)  The  vectors     -—^  e        'and  ■      e         J     attain  definite  limit  values 

A  (p,q)  and  3  (p,q)   if  infinity  is  approached  along  a  ray  p,   q,   r  of  the  instrument. 
The  vector  A   (p,q)   is   identical  with  the  vector  A  v/hich  has  been  introduced  in 
(2.2U)   in  order  to  characterize  the  geometric  optics  of  the  instrument.      In  other 
words   it    is   assumed  that   at   infinity  both  the  actual   electromagnetic   field  and  the 
geometric   optical  approximation  have  the  same  asymptotic  behavior. 

By  application  of  Green's   theorem  it   can  be   shovm  that   the   above  hypo- 
thetical solution  can  be   obtained  by  superposition  of  plane  v/aves,   namely  by  the 
integrals 


{f  =     JL  // ,^(p   q)    e^^^P^  +  qy  +   rz  +  W(p,q)       dpdq 

r  Vl-P     -  q  (2.32) 


->  1 

V   = 


2^//©^(p,q)    e^^^P-^  +  qy  +   rz  +  ¥(p.q)        dpdq 


f,  2  2 

p  VI  -  p    "  q 

the   domain  of   integration,  f~"^  is   determined  by  the   aperture  of   the   instrument.      It 

2  2 

IS  a  domain  of  the  p,    q  plane  which  lies   inside   the  unit   circle  p     +  q     =1..     The 
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integrals  (2.32)  axe  solutions  of  the  equations  (2.12)  for  any  choice  of  vectors 
A,  B  provided  thp.t  the  conditions 

IP  X  B  +  i.   =0 


->  -^    ->        ^  (2.33) 

]»  XA  -  B   =0;   p  =  (p.q.r) 


are  satisfied. 


2.U  The  evaluation  of  the  integrals  (2.32)  by  knovm  functions  is 
possible  only  in  a  fev/  sinple  cases.   Thus  \ie   r.re  led  to  the  problem  of  determin- 
ing the  principal  terms  of  these  integrals  by  an  asymptotic  development  with 
respect  to  a  parameter  proportional  to  the  quantity  k  =  r—  .   The  solution  of 

this  problem  then  should  provide  the  ansuer  to  the  question  whether  the  integrals 
(2,32)  solve  the  above  formulated  problem  of  correspondence  to  a  prescribed 
situation  of  geometrical  optics.   The  ansv/er  will  be  in  the  affirmative:   The 
principal  terus  of  the  integrals  (2. 32)  are  given  by  the  expression 


„  v/iT   ik  U/ 
u_  =  U  V —  e   r 


°      "  (2.1+1) 


A    ik^ 


V  =0V  ^^ 

0      k 

provided  that  the  dinensionless   quantity 

This  result  shoivs   that   an  asymptotic  development  beginning  v/ith  (2.UI)  must  break 
down  in  the  neighborhood  of  caustics   or  focal  points  where  K  =  00  . 

2,5     We  approach  the  problem  of  asymptotic  development  as  follows; 
ViTe  introduce  the  vectors 


?(x.y.z.t)   =  ^ff  ^(P'q^        dpdq 


V^    2    2 


9^<ct     yi-p  -q 


(2.51) 


a(x.y,3.t)    =^ff  ^^P'^^         dpdq 

</^<ct    v/l-p2-q2 


\ 
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(f^ct 


Fig.  h 
integrals  (2,32)  in  the  fonn  of  simple  integrals: 

-.00 

"o 
.00 


The  domain  of  integration  is  that 
part   of  the  domain  /where 

0''  =  px  +  qy  +  rz  +  W  <  ct 
We  then  may  write  the  multiple 


^  =  -  ic  /      i?«(t)  e 
00 


^    =-1= 


ot 


dt 


dt 


(2,52) 


The  vectors  P'  and  '^'  are,  at  any  fixed  point  x,y,z,  different  from  zero  only 

in  a  finite  interval  T,  <  t  <  T  ,  determined  ty  the  extreme  values  of  the  ftinction 

—  Cfxii   the  domain  /~^ ,   In  general  P  and  Q  v/ill  depend  analytically  on  t;   therg 

are  hov/ever  certain  points  t.^  inside  the  above  interval  where  P  and  Q,  or  their 
derivatives  are  discontinuous.  These  points  can  "be  formed  "by  studying  the 
topology  of  the  level  lines  ^f  =   const,  in  /^ , 

Consider  for  example  a  point  x,y,z  in  the  regular  part  of  the  field 
of  rays  (i.e.,  only  one  ray  passes  through  x,  y,  z).   In  this  case  ^possesses 
exactly  one  stationary  point,  P  ,  in  the  interior  of  /  ;  The  associated  value, 
ct  ,  is  either  a  maximum  or  minimum  of  6Cand  equa.ls  the  optical  path  <i/(x,y,z) 

from  the  source  to  the  point  x,y,z.  The  level  lines  (P  =   constant  thus  surround 

the  stationary  t»oint  P  in  the 
"      0 

manner  indicated  in  Fig  5.  There 
exists  two  level  lines 

f  =ct^  .  Lf=   ct^ 
(perhaps  more)  which  are  tangential 
to  the  boundary  of  /"*  at  points  P-,  and  P  .  We  can  find  these  level  lines  hy 

determining  the  stationary  points  of  the  boundary  values  of  the  function  (^  .  If 
A  and  B  are  analytic  in  /^then  ptt)  and  (^(t)  v/ill  he  analytic  functions  of  t 


except  at  the  points 


'0=  k'f'    \-     J^=   *2  =  |V^2 


1^. 


-^ 


given  by   the  points  P    ,   P    ,   P     of  Fig.    5«      '^^^  components   of  the  vectors  P     and 
Q,     thus   are   functions   of   t   of  the   type   sho\m  in  Pig.    6. 


to     ^.     ^X 

Fig.  6 

If  x,y,z  lies  outside  the  field  of  the  rays  then  there  is  no  interior 

stationary  point  of ^  but  again  two  or  more  stationary  boundary  points  P^  and  P  . 

p  Hence  there  are  only  two  points 


v/here  P  and  ^  are  irregular. 


Fig.  7 


If,  finally,  x,y,z  lies  in  the  irregular  part  of  the  ray  system,  i.e.,  if 
more  than  one  ray  passes  through  x,y,2;,  then  there  v/ill  be  more  than  one  station- 


ary point  in  the  interior  and  a  corresponding  number  of  stationary  boundary  points. 
In  this  case  there  exist  a  certain  finite  immber  of  points  tyu  ~  ^  luS'^t^t^)   where 

->   ->  r  r 

P  and  <^  are  non-analytic, 

— >    — > 

2.6  ^et  us  nov/  assume  that  P  sind  Q,  are  sectionally  analytic  functions 

of  t,  the  different  branches  being  separated  by  a  discrete  finite  set  of  points 

We  exclude  these  points  by  arbitrarily  small  intervals  and  transform  the  expressions 
(2.52)  by  repeated  partial  integration.  We  recognize  that  u  and  v'  admit  an 
asymptotic  development  of  the  type 

N 


t 


^ 


/    U^(x,y,z)  e 

/u=i  r 


ik 


^^ 


(2.61) 


N 


/>=1 


V^(x,y,z)  e 


ikl/^ 
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where  the  functions  "U u.   and  Vu,   depend  only  on  the  values  of  P  (t)  and  Q  (t)  in 
the  immediate  neighborhood  of  t  =  tu,  =  —   yytxCx^y,  z)  . 

The  asymptotic  series  of  IL^  and  Vuv/ith  regard  to  the  parameter,  k,  de- 
pends on  the  type  of  singularity  of  the  vectors  E  and  Q  at  the  points,  t  ^^  ;   It 

/ 
is  possible  to  derive  the  following  theorem: 

Let  us  assume  that  P*(t)  in  the  neighborhood  of  t  =  ti^ allows  a  representation 
of  the  following  form 

p'(t)    =  (t^-t)     YZ  ^n^V'**^''  ;  t  <  t^ 

.      ^  (2.62) 

p'(t)  =  (t  -g'    Y2  V^-y^""  ;  t  >  t^ 

Then  the   function  ILj  in  (2.6l)  may  be  developed  into  an  asjinptotic  series  of  the 
fonn 

CD        q^/~^(n+l+/3)  P^  P(n+1+X)1 


V  =  a 


2.7  ^^s  hp.ve   reduced  the  problem  of  asymptotic  evaluation  of  the 
diffraction  integrals  to  the  problem  of  finding  the  character  of  the  irregulari- 
ties  of  the  vectors  P  and  Q,  at  certain  singular  points  t^.  These  singular 
points  are  determined  by  the  system  of  the  level  lines  (J)    =   constant  In  f~^  ,   in 
particular  by  the  stationary  points  of  the  interior  or  of  the  boundary  values  of 
the  phase  (J),     Me   mention  that  an  interior  regular  maximum  or  minimum  point  of  (/? 
introduces  an  ordinary  discontinuity  {ff^=lS=   0)  of  p'^.t),  'a   stationary  boundary 
point,  hov/ever,  a  square  root  discontinuity  <X=  I/2  or  P  =  I/2,  This  shows  that 
the  asymptotic  contribution  associated  with  stationary  boundary  points  (P.  and  P 

in  Fig.  5)  a^e  of  higher  order  in  —  than  those  from  interior  stationary  points. 

One  also  finds  that  the  first  coefficient  P  or  Q  .in  case  0<  =  P=  0,   has  the 

00 

value  k  y/K^     so  that  the  principal  term   in  the  development   (2,6l)  has   the  form 

°    ^  (2.71) 


as  was   stated  above  in   (2.UI), 


0  k  '' 
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If  the  function  ^ possesses  more  than  one  interior  stationary  point  then  a 
finite  sum  of  terms  (2.71)  is  obtained  as  principal  asymptotic  contritiution.  This 
is  the  case  if  more  than  one  ray  passes  through  the  point  x,y,z  in  question. 

Now  it  can  he  shown  that  the  veciors  P  and  (^  are  indeed  the  pulse  solutions 
of  Maxwell's  Equations  for  this  particular  prohlem.  The  general  theory  requires 
that  one  kno\/  the  pulse  solution  in  order  to  determine  the  asymptotic  hehavior  of 
the  steady  state  field  and  this  ptilse  solution  is  known  explicitly  in  the  form  of 
the  vectors  P  and  Q  for  the  particular  protlera  of  transmission  through  a  lens  that 
has  heen  presented  in  this  sketch, 

3.      There  are  several  other  items  of  research  planned  under  this  general 
study  of  asymptotic  hehavior  of  steady  state  fields.   It  may  he  possible  to  show 
Lhat  the  general  theory  described  under  the  first  item  above,  namely  the  relation- 
ship between  pulse  solutions  and  steady  state  fields  can  be  approached  through 
general  integral  relationships  v/hich  v;ill  avoid  the  difficulties  now  associated 
with  the  intricate  use  of  complex  variable  theory.  A  second  possible  application 
and  one  long  envisioned  is  to  the  problem  of  anomalous  propagation,  where  the  cal- 
culation 01  the  fields  can  be  represented  by  integrals  of  the  same  form  as  those 
occurring  in  the  theory  described  above, 

V,  Inversion  I-iethod  and  Perturbation  Method 

1.      The  problem  of  propagation  in  a  spherically  symmetric  medium  was  con- 
sidered in  an  early  Report  \2]   for  a  very  special  variation  of  index  of  refraction, 
It  was  assiuned  that  the  atmosphere  of  the  earth  consisted  in  a  layer  of  finite 
height  in  \rhich   the  index  of  refraction  varies  inversely  as  the  square  of  the  dis- 
tance from  the  center  of  the  earth.  It  is  assui.ied  that  above  this  layer  the  index 
is  constant.   The  earth  is  supposed  poorly  conducting. 

It  v;as  recognised  that  the  choice  made  of  refractive  index  in  the 
atmosphere  did  not  correspond  to  physical  reality  but  it  was  hoped  that  a  method 
similar  to  that  in  the  paper  might  be  applied  to  a  sufficient  variety  of  refractive 
indexes  to  yield  useful  conclusions  for  the  problem  of  anomalous  propagation.  This 
hope  was  not  borne  out  and,  in  spite  of  the  mathematical  interest  associated  \7ith 
certain  aspects  of  the  problem  approached  in  this  v;ay,  the  general  program  came  to 
be  abandoned.   In  an  early  Progress  Report  submitted  to  the  Laborato.:ies  there  is 
some  discussion  of  the  types  of  variation  in  refractive  index  which  can  be  handled 
by  the  method  of  the  paper,  and  some  numerical  calculations  shouixg  their  inadequacy 
to  the  problem. 


22. 


The  Report  {^z']  uses  the  classical  transformation  of  space  kno\m  as  in- 
version,  given  by  /^  =  a  /r,  where  r  is  the  distance  of  a  point  from  the  origin, yO 
the  distance  to  the  origin  of  the  transformed  point  collinear  v/ith  the  origin  and 
the  first  point.  The  constant  a  represents  the  radius  of  the  inversion-sphere,  in 
this  case  the  radius  of  the  earth.  The  vave   equation  satisfied  "by  the  Herzian 
potential  function  in  the  case  of  constant  index  of  refraction  is  transformed,  "by 
inversion,  to  an  equation  corresponding  to  an  index  varying  according  to  the  in- 
verse square  law.  The  complete  solution  of  the  transformed  equation  may  "be  derived 
from  the  known  solutions  to  the  original  equation.  There  remains,  however,  the 
difficulty  of  finding  solutions  consistent  with  the  proper  boundary  conditions.  A 
solution  to  the  problem  is  obtained  in  the  form  of  a  slowly  convergent  series. 

2.      The  problem  of  propagation  in  non~homogeneous  media  was  investigated, 
for  a  time,  from  the  point  of  view  of  certain  modifications  of  classical  per- 
turbation methods.  This  attempt  is  discussed  in  a  preceding  Qjiarterly  Progress 
Report, 

VI.  Wiener-Hopf  Methods  in  Propagation  Problems. 

In  past  years,  considerable  progress  had  been  made  by  members  of  the 
Radiation  Laboratory,  located  at  a.  I.  T.,  in  the  aj^lication  of  integral  equations 
of  the  Wiener-Hopf  type  to  problems  in  wave-guide  and  in  propagation  theory.  It 
has  not  always  been  easy  to  know  v/hat  class  of  problems  remained  that  might  be  ex- 
ploited by  this  or  related  methods.  In  connection  with  a  recent  expansion  in  the 

forfner 
program  undertaken  by  the  group,  one  of  the  active  »embers  of  the  Jjladiation  Labora- 
tory v;as  engaged  as  consultant,  and  this  class  of  problems  is  now  being  investigated 
by  members  of  the  group.   So  far,  the  activity  has  been  primarily  that  of  study  and 
group  discussion  in  a  seminar  devoted  to  this  topic, 

VII.  Bibliography  ' 

In  the  course  of  study  and  exchange  of  ideas  by  members  of  the  group, 
and  through  contact  with  other  groups  engaged  in  related  research,  there  began  to  be 
collected  a  substantial  Bibliography  on  the  general  topic  of  the  propagation  and 
scattering  of  electromagnetic  waves.  This  was  subsequently  augmented  by  a  systematic 
search  through  periodicals  and  references  appearing  in  recent  years.   It  is  sub- 
mitted as  a  Report  [9j  . 


23. 

VIII.  Education 

Serious  attention  was  given  to  one  of  the  objectives  of  government  sub- 
sidization of  university  research,  namely,  the  development  of  the  scientific  poten- 
tial of  the  country.   Considerable  training  in  the  field  of  electromagnetic  radiation 
v/as  afforded  a  few  yoiinger,  pre-doctorate  mathematicians,  and  much  interest  in  the 
application  of  mathematics  to  electromagnet ism  stimulated  among  the  students  and 
faculty  at  New  York  University. 

Two  offerings  in  the  Graduate  School  also  aided  in  stimulating  interest 
in  the  applications  of  mathematics  to  electromagnetic  theory.   During  the  academic 
year  19^+6-^7  a  seminar  in  electromagnetic  problems  and  during  the  year  19147-1+8  an 
advanced  course  in  electromagnetic  theory  were  offered.  Both  v/ere  v;ell  attended. 
While  these  courses  \;ere  not  subsidized  by  government  funds  the  incentive  to  offer 
them  and  the  interest  in  them  on  the  part  of  some  of  the  students  was  motivated  "by 
association  with  the  work  of  this  group, 
IX. _ Administrative  Report  (for  1  July  to  3O  September) 

1.  Personnel 

Dr.  Nathan  Marcuvitz  has  been  appointed  to  this  project  as  consultant. 
His  services  are  on  a  part  time  basis,  commencing  1  July,  I9I+8.  Dr.  Marcuvitz  is 
Assistant  Professor  at  Brooklyn  Polytechnic  Institute,  and  has  a  background  of 
active  participation  on  the  theoretical  staff  of  the  Reuiiation  Laboratory. 

Dr.  Samuel  Karp  has  been  appointed  to  this  project  as  senior  scientist* 
His  services  are  on  a  fhll  time  basis,  conirnencing  1  August^  I9U8,  Dr.  Karp  has  had 
extensive  research  experience  at  Brown  ITniversity  in  mathematical  techniques  which 
we  are  seeking  to  apply. 

2,  Correspondence 

Security.  On  7  July,  19l|g,  we  received  memo  3U5-IOOO-2  regarding 
security  questionnaires  from  the  Air  Materiel  Command,  which  \/e  acknowledged  re- 
questing a  supply  of  forms.  On  I3  July  Lt.  Col.  M.C.  Edenfeld  called  our  attention 
to  paragraph  2  subparagraph  (b)  of  the  Annual  Secrecy  Agreement  and  on  11  August 
Col.  Theodore  U.  Hatt,  of  Watson  Laboratories,  confirmed  our  understanding  that 
in  the  case  of  non-aliens  seciirity  clearance  is  required  only  for  participation  in 
projects  of  "secret"  or  "top-secret"  classification. 

On  23  August  ve  \/ere  advised  by  Col.  Theodore  M.  ilatt  that 
clearance  is  requested  for  subcontractors,  where  required  under  article  20. 
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We  replied  to  lir.  S.S.  Hosenterg,  of  V/atson  Laboratories,  26  August,  that  Doris  A. 
Troy,  ^-5  Astor  Place,  ilew  York,  and  also  the  Duplicating  Department  of  Hew  York 
University  do  nimeographing  for  us.  Correspondence  is  in  progress  to  ascertain 
v;hether  this  constitutes  subcontracting,  in  the  meaning  of  the  article  20, 
Property 

On  1  September  ve   wrote  to  Lt.  Henry  W.  Liljedahl,  referring  to 
earlier  coi.inunications,  in  connection  with  the  transfer  of  property  charged  to 
Contract  Ko.  W28-099-a-c~170,  and  received  reply  that  this  property  has  been  trans- 
ferred to  this  Contract,  ITo.  './28-099--ac-17 2, 
Miscellaneous 

On  5  July  i;e  informed  lir.   xT.  C,  Gerson  that  Dr.  '/ilhelm  Magnus, 
whose  services  for  this  project  we  had  earlier  solicited,  was  being  brought  to  this 
country  in  connection  i/ith  certain  Navy  contracts. 

On  20  July  v;e  v;ere  advised  by  the  Air  Materiel  Command  that  the 
University  of  Michigan  had  requested  tv/o  copies  of  Research  Report  iTo.  172-7.  en- 
titled "Reflection  and  Transmission  of  Electromagnetic  './aves  by  a  Spherical  Shell", 
by  Herbert  3.  Keller  and  Joseph  B.  Keller,  Permission  was  also  requested  to 
index  the  title,  author,  and  abstract  in  the  Air  Technical  Index  (A.T.I.)  catalog 
cards.  Me   forv/arded  a  copy  of  the  above  report  to  the  Air  Materiel  Command  on  3O 
July  and  expressed  our  desire  to  have  this  document  indexed. 

Me   have  received  Progress  Reports  of  the  Canputation  Laboratory  of 
the  national  Applied  Mathematics  Laboratories  for  the  months  of  June,  July  and 

August.   Reports  are  also  beirig  received  fro^.i  the  Cruft  Laboratory  of  Harvard 
University  and  the  Research  Laboratory  of  Electronics  at  i'i.I,T. 
3,     Conferences 

Several  members  of  the  group  attended  the  Second  Annual  Colloquium 
on  Applied  Mathematics,  held  in  Cambridge  at  li.  I.T,,  29-31  July.  The  Colloquitun 
this  year  was  devoted  entirely  to  electromagnetics.  Dr.  R.  laineberg,  of  our  group, 
was  one  of  the  invited  speakers,  unfortunately  not  able  to  participate  by  reason  of 
illness. 

Dr,  Paul  M.  Marcus  of  the  Office  of  ITaval  Research,  visited  this 
group  on  10  August,   We  discussed  the  general  nature  of  our  v/ork.   It  is  hoped  that 
Dr.  Marcus  v/ho  will  travel  in  England  on  behalf  of  the  O.IT.R.  may  be  able  to  assist 
us  in  establishing  an  interchange  of  ideas  v/ith  v/orkers  in  England  on  projects 
similar  to  ours. 
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Dr.  Marvin  Chodorov  of  Stanford  University  visited  us   on  I5  August. 
We  discussed  the   relationship   of  our  \/ork   to   the   work  teing  done  at   Stanford 
University  under  the   direction  of  Professor  W.   'i.   Hansen, 

X.   Research  Reports  Uos.    172~1  to  172-9    sutoitted  under  Contract  ITo.W'2g-099-ac-172 

1.  Reflection  of  Electromagnetic  Waves,   "by  Joseph  B.   Keller. 

2.  An  Application  of  Inversion  to  \<'ave  Propagation,   fey  Bernard  Friedman. 

3.  The  Propagation  of  Slectroaagnetic  Plane  Waves   in  Plane  Parallel  Layers, 
"by  R.   K,   Luneterg. 

U.  Propagation  of  Dipol3  Radiation  Through  Plane  parallel  Layers,   by  Jerome  Lurye. 

5.  Reflection  and  Transmission  of  Slectromagnetic  Waves  Tjy  Thin  Curved  Shells, 
Joseph  B.   Keller. 

6.  Propagation  of  Electromagnetic  V/aves  from  an  Arbitrary  Source  Through  In- 
homogeneous  Stratified  Atmospheres,   by  R.   K,-  Luneberg. 

7.  Reflection  and  Transmission  of  Electromagnetic  V/aves  by  a  Spherical  Shell, 
by  Herbert  3.  Keller  and  Joseph  B.   Keller, 

8.  Maxwell's  Equations   in  Spherically  Symmetii.c  Media,   by  R,   K.   Luneberg. 

9.  Bibliography:     Propagation  and  Scattering  of  Electromagnetic  Waves,   compiled 
by  Thelma  Braverman  and  Jerome  Lurye. 
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